A supermassive binary black hole with triple disks 
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ABSTRACT 

Hierarchical structure formation inevitably leads to the formation of super- 
massive binary black holes (BBHs) with a sub-parsec separation in galactic nuclei. 
However, to date there has been no unambiguous detection of such systems. In 
an effort to search for potential observational signatures of supermassive BBHs, 
we performed high-resolution smoothed particle hydrodynamics (SPH) simula- 
tions of two black holes in a binary of moderate eccentricity surrounded by a 
circumbinary disk. Building on our previous work, which has shown that gas 
can periodically transfer from the circumbinary disk to the black holes when the 
binary is on an eccentric orbit, the current set of simulations focuses on the for- 
mation of the individual accretion disks, their evolution and mutual interaction, 
and the predicted radiative signature. The variation in mass transfer with or- 
bital phase from the circumbinary disk induces periodic variations in the light 
curve of the two accretion disks at ultraviolet wavelengths, but not in the opti- 
cal or near-infrared. Searches for this signal offer a promising method to detect 
supermassive BBHs. 

Subject headings: black hole physics - accretion, accretion disks - binaries:general 
- galaxies:nuclei 
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Introduction 



1.1. Scientific Motivation 



Supermassive black holes are common in galactic nuclei and are considered to have coe- 
volved with their host galax i es (Magorrian et alii 998 ; Ferrarese &: Merritt 2000 ; Gebhardt et al 
2000l ; lYu fc Tremaind 120021 ; |Pi Matteo et al.ll2005l ). Since galaxies gr ow through mergers, i t 
is inevitable that supermassive binary black holes (BBHs) will form (jEbisuzaki et al.lll99ll ). 
Supermassive BBHs evolves via three stages (IBegelman et al.lll980l ; lYull2002l ). First, each of 
the black holes sinks independently toward the center of the common gravitational potential 
due to dynamical friction with neighboring stars. When the separation between the two 
black holes becomes less than ~ 1 pc, angular momentum loss by dynamical friction s lows 
down due to the loss-cone effect (jRooslll98ll ; lMakinolll997l ; iMilosavljevic fc Merrittll200ll ). A 
hard binary forms during the second stage, whose orbit slowly shrinks, until the final stage 
where the black holes merge by emitt ing gravitational ra diation. Recent hydrodynamic sim- 
ulations of a gas- rich galaxy merger (IMayer et al.l 120071 ) show that the interaction between 
the black holes and the surrounding stars and gas rapidly leads to the formation of a BBH. 
After the two black holes have reached the scale of the circumnuclear disk (~ 100 pc), it 
takes less than 1 Myr for them to form a hard binary on scales of ~ 1 pc. 

Despite the general expectation that BBHs should form naturally during the course 
of galaxy formation and evolution, the currently available observational evidence for BBHs 
have been largely circumst antial and indirect. T hese include the detection of X-shaped 
structures in radio galaxies (IMerritt fc Ekersll2002l ). double compact cores with a flat radio 
spectrum (Rodriguez et al. 20061 ) . and close pairs of AGNs such as that seen in NGC 6240 
( IKomossa et al.l 120031 ) . As these systems evolve, we expect that their black holes would 
eventually form a strongly gravitationally bound binary on sub-parsec scales. Such small- 
scale bin aries have been inferr ed to be present from the quasi-periodic outbursts detected in 
OJ 287 fail 



anpaa et al. 



in 3C 66B (ISudou et al. 



1988) and the proper motions seen toward the compact radio core 



20031 ). but the interpretation of neither case is unambiguous. 



In view of the astrophysic al importance of BBHs , it would be desirable to devise robust 
methods for detecting them. ISundelius et al.l (119971 ) showed that when two unequal-mass 
black holes interact, for instance in the aftermath of a minor galaxy merger, the impact of 
the secondary black hole on the accretion disk around the primary black hole may generate 
periodic flares. These authors invoked this model to explain the brightness variations in 
OJ 287. In supermassive BBHs with extreme mass ratios, the secondary black hole could be 
embedded in the accretion disk of the primary black hole and migrate toward it. The kind 



of disk-binary interaction offers a promising solution to the loss-cone problem (jlvanov et al. 
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19991 : iGould fc RbJhoOol : Urmitage fc Nataraianll2002l . l2005h . H (|200lh proposed that BBHs 
might be detectable from the integrated Fe Ka line of two accretion disks, which would give 
rise to an unusual profile with double or multiple peaks. 



1.2. Our Previous Work 



It is likely that a rotating gas disk (i.e., circumbinary disk) forms around supermassive 
BBH during the second stage of a gas- rich galaxy merger. If so, what will happen as a r esult of 
the interaction between the BBH and the circumb inary disk? Havasaki et al.l (120071). using 



a sm oothed particle hydrodynamics (SPH) code (IBenz et al. 



199C 



u~ BenzJ 



1990 



Bate et al. 



19951 ). investigated the accretion flow from a circumbinary disk onto the supermassive BBH 
on sub-parsec scales. In that study, the circumbinary disk was aligned with the binary 
orbital plane. They found that mass transfer can occur from the circumbinary disk to the 
supermassive BBH in two steps. First, the initially circularized circumbinary disk becomes 
elongated due to the azimuthal m = 2 component of the binary potential. The gas density 
then grows at the two points closest from the black holes. Next, when the gas reaches beyond 
the potential barrier of the binary, mass inflow is initiated by gas pressure so that the gas 
freely inspirals onto each of the black holes. This process repeats every binary orbit. The 
viscosity of the circumbinary disk is ineffective for the mass transfer process because the 
viscous timescale is much longer than the orbital period. For BBH systems on sub-parsec 
scales, the viscous timescale of the circumbinary disk can be written as 
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where ass = 0.1 is the Shakura & Sunyaev (1973) viscosity parameter, Mbh is the total black 
hole mass, and a is the semi-major axis, respectively. The inner edge of the circumbinary 
disk is assumed to be at the 1 : 3 resonanc e (~ 2.08a), where the visco us force is balanced 
with the tidal/resonant force of the binary (jArtymowicz &: Lubowlll994l ). The circumbinary 
disk is assumed to follow an isothermal equation of state. 

The mass transfer rate significantly depends on the binary orbital phase in eccentric 
binaries, whereas it shows little variation with orbital phase in circular binaries. In both cases 
the supermassive BBH system is expected to have triple disks because the mass transfer from 
the circumbinary disk would inevitably lead to the formation of an accretion disk around 
each black hole. However, there is little known about how the accretion disks are formed 
and what signatures they exhibit. We have, therefore, performed a new set of simulations 
at higher resolution, adopting the same set of binary orbital parameters we had previously 
used except for a semi-mjor axis a = 0.01 pc, which was chosen to yield a shorter period 
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binary than that of the previous simulation (a = 0.1 pc; Hayasaki et al. 2007). 

The simulations should ideally take into account all of the processes at work in a triple- 
disk system, including mass transfer from the circumbinary disk to the individual black holes, 
and the accretion flow onto each black hole. Such a simulation, however, would require an 
enormous amount of computational time. Therefore, we divide these processes into a two- 
stage simulation. In the first stage, we investigate how the mass is transferred from the 
circumbinary disk to the effective gravitational radius of each black hole. In the second 
stage, we confine ourselves to simulate only the accretion onto each black hole by setting 
the mass-transfer rate from the first-stage simulation as a new boundary condition. In this 
paper, we focus on the results of the second-stage simulation. 



2. Numerical Procedures 

We set the binary orbit on the x — y plane with the semi-major axis along the x-axis. 
Initially, each black hole is at the periastron. It orbits with P or b ~ 9.4 yr, and the orbital 
eccentricity is chosen to be e = 0.5. As a result of the interaction between the black holes 



and the neighboring stars and gas (jMayer et al.l 120071 ). the BBH likely forms rapidly on 



scales of a few parsecs. The total black hole mass is 1O 8 M , which is a typical value for the 
centers of massive elliptical galaxies. The simulations were carried out with a total of 152,000 
particles at the end, each one having a mass of 1.0 x 10~ 15 Mbh- For simplicity, we set the 
mass ratio to q = M2/M1 = 1.0 and assume that the circumbinary disk is coplanar with the 
binary orbital plane. The effect of the mass ratio, the orbital eccentricity, and inclination 
angle on the mass transfer process will be studied in a subsequent paper. In this section, we 
first explain the boundary conditions given by the first-stage simulation, followed by a brief 
description of the energy equation. 



2.1. Boundary Conditions 



To reduce computational time and to significantly improve the resolution, we confine 
the simulations only to the accretion flows around the supermassive BBH by supplying mass 
periodically from the outer boundary conditions, which are defined by the SPH particles 
captured by the black holes (jHayasaki fc Okazaki 2004 ; Hayasaki et al. 2007 ). The capture 
radii correspond to the effective gravitational radii 0.8 tl, where ri is the innermost common 
gravitational radius of the binary potential for a circular binary. The mass is added at a rate 
of 1 M Q yr -1 on the outer edge of the circumbinary disk. The system reaches quasi-steady 
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state equilibrium by 30 binary orbits, when the mass input balances the mass output. The 
mass of the circumbinary disk is ~ 10~ 6 Mbh, and there are 56, 870 SPH particles at the end 
of the simulation. To reduce the fluctuation noise, the data are folded over the orbital period 
30 < t < 60, where the unit of time, unless otherwise noted, is P or b — 9.4 yr. The orbital 
phase dependence of the mass transfer rate obtained by this procedure is shown by the solid 
line in Figure 1. We model the mass transfer process from the circumbinary disk to the 
supermassive BBH by injecting gas particles at a given phase- dependent rate with spatial 
and velocity Gaussian distributions derived from the averaged values over 30 < t < 60. 



2.2. Energy Equation 

The new simulations no longer assume an isothermal equation of state. We introduce 
an energy equation in order to resolve the temperature profile. The specific internal energy 
of an ideal gas is explicitly given by 

u = -^TP 7 -\ (2) 
7 - 1 

where p is the density, 7 is the adiabatic index for an ideal gas (5/3), and K is a function of 
the entropy, which is obtained by the first law of thermodynamics. Time differentiation of 
equation ([2]) yields the energy equation, which can be written as 

du P dp ^ p 7-1 dK 
dt p 2 dt 7 — 1 dt 

The first term on the right-hand side gives the heating of the gas at pressure P. In the 
current simulation, we assume that the heating due to viscosity is dissipated by cooling due 
to blackbody radiation. The cooling term is related to the effective temperature T by the 
Stefan-Boltzmann law, 

2aT 4 7 -i dR 
pH c 7 — 1 at 

where a is the Stefan-Boltzmann constant and H c is the characteristic scale of the system, 
and the factor 2 represents radiation from the two sides of the accretion disk. 

In SPH calculations, gas heating is implemented using a standard SPH artificial viscosity 
jMonaghan fc GingoldlflQsl ) 



ij \ for ■r ij >0 U 

for particle i and its neighboring particle j, where osph and /9sph are the linear and nonlinear 
artificial viscosity parameters, respectively, p.y = (pi+pj)/2 is the mean density, 7*y = r^ — Tj 
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is the relative position vector, Vi 
hijVij ■ r y /(r y + rjij) with i/'j 



0.01/i?- and 



Vj is the relative velocity vector, and fiij = 
(hi + h j) 1 2 being the mean smoothing 



length. The entropy change of particle % can be written as (jBatd 119951 ) 



dKi 



7 



2p] 



-y N n 

— 22 "'; 11 ';'' 7 './ ' ^iW^Vij, hij), 



(6) 

where N n is the number of neighboring particles and W, proportional to \hij\~ 3 , is called 
the SPH kernel. Thus, the effective temperature of particle i can be calculated from equa- 
tions (SD-dSD: 



Ti 



1/4 



(7) 



where we set H r 



h^ in the current simulation. 



3. Accretion Disks around the Black Holes 

Each of the two black holes has an accretion radius 5 x 10~ 3 a, which is the radius of 
the inner boundary of the simulation and corresponds to r m ~ 10.4 tbh, where tbh is the 
Schwarzschild radius of a black hole with mass 5.0 x 10 7 M Q . To verify the robustness of our 
simulation results, we have compared the averaged mass accretion rate during one orbital 
period with that of a simulation whereby the number of injected particles was doubled and 
the particle mass was reduced in half. We find no qualitative or quantitative differences to 
within ~ 8%. 



3.1. Formation Process 

We use the standard SPH artificial viscosity throughout the simulations. In simulat- 
ing the disks, the shear component of the SPH artificial viscosity can be replaced with 
the Shakura-Sunyaev viscosity parameter ass i n the following form: ass = O.laspH^-/ H 



(IHayasaki et al.l 120071 ). where h is the azimuthally averaged smoothing length and H is the 
azimuthally averaged scale height of the disk. The estimated value of ass ranges from 0.001 
to 0.2 over the entire region of the disk. Little gas can accrete viscously onto the black holes 
within the simulation run time because the viscous timescale is much longer than the orbital 
period (cf. Eq. |SJ). 

The injected particles originally have eccentric orbits around the black holes because 
their angular momenta are smaller than that of the particles at the disk outer edge (0.2 a). 
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Such gas particles accrete directly onto the black holes within one dynamical timescale of the 
disk, T dyn = 27r/f2 disk ~ 0.1P or b, which is estimated at the disk outer edge. In addition, the 
particles in each disk lose angular momentum mainly from tidal interaction with the other 
black hole because the timescale for angular momentum loss at the outer edge of the disk is 
shorter than the binary orbital period: r t id c ~ 27r/(f2disk ~ ^orb) ~ O.lPorb; where f2di S k and 
f2 or b are the angular frequency of the disk and of the binary, respectively. This also makes it 
possible for the gas to accrete onto the black holes within the timescale of the binary orbit. 

We illustrate the behavior of the two accretion disks in Figure 2, which gives snapshots 
of the accretion flow around the supermassive BBH covering one entire orbital period during 
the 10th binary orbit. At the end of the simulation, the disk has a mass of ~ 10~ 8 Mrh- The 
two white crosses on each density map indicate the mass supply points. When the two black 
holes are closest together (t = 10.0), the disks are deformed due to the tidal interaction, 
and the supplied gas is added at the outer parts of the disks (t = 10.04). At t = 10.1, the 
gas particles that are added at the disk outer edge start to inspiral onto the black holes 
within a dynamical time. The disk size become ~ 50% smaller than the disk outer edge. As 
the black holes recede toward apastron (t = 10.18), the disks are restored from their tidal 
deformation and gradually circularize; tidal tails connect to the mass supply points between 
the two black holes. As can be seen in Figure 3 a, during this time the temperature of the 
inner part of the disk is highest. The disks are connected to each other via tidal tails and 
are still eccentric, even at apastron (t = 10.5). The disks approach each other again toward 
the next periastron, and the injected particles once again start to be added to the disks 
(t = 10.75). This pattern repeats periodically. 



3.2. Survival 

Whether the accretion disks formed around the black holes are transient or persistent 
depends on the viscous timescale of the disk compared with the orbital period. For simplicity, 
we assume the accretion disk to be geometrically thin and isothermal with a Shakura-Sunyaev 
viscosity parameter «ss- The ratio of the viscous timescale of a disk of size R& to the orbital 
period is then given by 



^ ~ 5, x !«,. ' m 1/2 m M » ■ (S> 



It is immediately seen that the viscous timescale is much longer than the binary period. 
Thus, the accretion disks will survive, not only over the whole orbital phase but also over 
the simulation run-time (~ 10 P or b) once they are formed. 
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If the separation between the two black holes is less than 0.01 pc, emission of gravi- 
tational waves becom es the dominant mechanism to dissipate the binary's orbital energy 
( IBegelman et al.lll980l ). Dissipation makes the binary evolve rapidly toward the coalescence 
of the two black holes. Therefore, we should consider whether the accretion disks survive 
with shrinkage of the semi-major axis due to gravitational wave emission. The merging 
timescale bet ween two black holes due to the emission of gravitational radiation is given by 



fjPetersll!964 ) 
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(9) 



where /(e) = (1 - e 2 ) 7 / 2 /(l + 73e 2 /24 + 37e 4 /96). From equations (7) and (8), we can 
compare the viscous timescale of the disk with the merging timescale: 

i/ 2 /ni\ /m4zr\ /n ni t^^\ 7 / 2 / M„ TT \ 7 / 2 

l>M . (10) 
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The range of our simulation parameters satisfy r gr > ^> P or b- This ensures that the 

accretion disks, once formed, would survive until the black holes coalescence if the mass 
supply from the circumbinary disk to the central binary can be maintained. Gravitational 
wave emission has little influence on the formation and evolution of the accretion disks. 



3.3. Temperature Profiles, Spectral Energy Distributions and Light Curves 

The radial distribution of the azimuthally averaged temperature is shown in Figure 3 a. 
The solid and dotted lines denote, respectively, the highest and lowest temperatures at r in . 
The gas is distributed over the entire region of the disk, and the temperature profile varies 
with orbital phase. Figure 36 shows the corresponding spectral energy distributions (SEDs) 
normalized by the Eddington luminosity for a total black hole mass Mbh = 1-0 x 10 8 M Q . 
The SEDs are computed numerically by integrating the spectrum, which is a function of 
disk temperature, from the outer edge of the disk (0.2 a) to r in . Each SED has a peak in the 
ultraviolet (UV) region. 

Figure 4a shows the orbital phase dependence of the light curves emitted from the disk 
in the UV, optical, and near-infrared (NIR) bands, obtained by averaging out the time vari- 
ations over the orbital periods 10 < t < 11. While the UV light curve shows substantial 
fluctuations with time, the optical and NIR light curves exhibit little variation. This is 
because a small temperature change produces large-amplitude variations in the UV region, 
where the SED decays exponentially (Fig. 36), whereas it has little effect on the Rayleigh- 
Jeans part of the spectrum (optical and NIR). What does the X-ray light curve look like? 
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The dashed line in Figure 1 shows the orbital variation of the mass accretion rate, which 
is defined as the number of SPH particles per unit time captured at r in over the orbital 
period 10 < t < 11. If the gas in the disk accretes onto the black hole while keeping this 
orbital phase dependence, the X-ray light variation is expected to track the variation in mass 
accretion rate. This situation seems possible, if the accretion timescale in the inner region is 



shorter than the binary orbital period 
radiatively inefficient accretion flo ws ( 



i.e., close to the free-fall timescale), as is the case for 



Kato et al.lll998l). T he recent hydrodynamical sim- 



ulations of supermassive BBHs by iBogdanovic et al.l (120081 ) predict strong X-ray outbursts 
during periastron passage. However, since our simulations do not take into account the rel- 
ativistic correction of the gravitational potential in our simulations, the detail properties of 
the X-ray emission are beyond the scope of this work. 

In order to isolate the effect of the mass transfer from the circumbinary disk onto the 
accretion disks, we carried out a simulation in which the mass transfer was artificially stopped 
over the orbital period 11 < t < 12. For comparison purpose, the light curves for this case 
are shown in Figure 4b. Note that none of the light curves varies much with orbital phase. 
Thus, the periodic mass transfer from the circumbinary disk is clearly responsible for the 
significant variation of the UV light curve. 



4. Summary and Discussion 

We have performed SPH simulations of accretion flows around each black hole in a 
supermassive BBH system, taking into account periodic mass transfer from the circumbinary 
disk. An accretion disk forms and survives around each black hole. Thus, supermassive 
BBHs have triple disks, consisting of two accretion disks, one around each black hole, and a 
circumbinary disk surrounding them as a mass reservoir. In such systems, three dynamical 
processes are expected to give rise to structures in the two accretion disks. First, periodic 
mass transfer makes the disks eccentric because the gas particles from the circumbinary disk 
originally have elliptical orbits around the black holes. This makes it possible to accrete 
directly onto the black holes during the binary orbit. Second, once the disks form, periodic 
mass transfer from the circumbinary disk episodically gives an impact to the outer edge of 
the disks, resulting in the excitation of a one-armed spiral wave in the disks (e.g., Hayasaki 
& Okazaki 2005). Mass exchange between the two disks further reinforces the one-armed 
spiral wave. Finally, the disks are tidally deformed by the time-varying binary potential 
owing to the orbital eccentricity. 

During the simulation, the accretion disks are still developing because their viscous 
timescale is much longer than the simulation run time (10P or b). In addition, a one-armed 
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spiral wave could also contribute to the non-axisymmetric disk structure because its propa- 
gation ti mescale, r wa , VP ~ (a^/27r)T^ dlsk ~ 10 3 P or b, is also much longer than the simulation 



run time (lHayasaki Sz Okazakill2005l ). In such an early stage of the disk evolution, the pe- 
riodic mass transfer has the dominant effect on the disk structure. Without mass transfer, 
tidal deformation alone has little influence on the light curve variation. As a consequence 
of the periodic mass transfer, we show that the light curve in the UV — and mostly likely 
also in the X-rays — exhibits variations with orbital phase, whereas little fluctuation is seen 
in the optical and NIR light curves. NIR variability might be detectable if the binary is 
surrounded by a dusty torus, whose inner edge would be irradiated by, and thus mirror, the 
variable central UV and X-ray sources. The photometric variability discussed in this paper 
is distinctive from the intrinsic variability normally seen from single accretion disks or jets, 
and thus provides a unique signature to identify BBH systems in galactic nuclei. 

The orbital period of the fiducial case simulated here, ~ 9.4 yr, is long for any realistic 
observational experiment. The only known BBH c andidate with a similar orbital period (11- 



12 yr) is the blazar OJ287 (jSillanpaa et al.lll988l ). Supermassive BBHs with shorter orbital 



periods (e.g., up to a few years) would be more practical observationally. However, the 
generic features of the triple-disk BBH system we have simulated would remain unchanged 
for shorter orbital periods. To avoid general relativistic corrections to the potential, the size 
of the accretion disk needs to be larger than ~ 30 tbh- Since the semi-major axis of an 
equal-mass binary is ^ 20 times the disk size, the results of our simulations can be scaled 
to systems with a ^ 600 tbh- This condition gives a wide range of tractable orbital periods, 
from P orb = 5 d for M BH = 10 6 M to P olb = 1.4 yr for M BH = 10 8 M . 

In the first-stage simulation, the orbital-phase dependence of mass transfer rate is 
smoothed out as the orbital eccentricity gets toward zero. On the other hand, the pea k 



of the mass transfer rate in a higher eccentricity binary is sharper (cf. lHayasaki et al.ll2007l ). 
These features can also be seen in the mass accretion rates and the light curves in the second- 
stage simulation. Most of the angular momentum gained from the mass transfer from the 
circumbinary disk is converted into that of the BBH via the tidal interaction between the 
BBH and the accretion disk. These processes damp the orbital eccentricity of the binary. 
However, since the damping timescale approximately equals the viscous timescale, which is 
much longer than the binary orbital period (cf. equation |Sj), this effect cannot be seen in 
this simulation. 

When the separation between two black holes falls below ~ 0.01 pc, the binary wi ll suffer 



rapid evolution due to the emission of gravitational radiation (IBegelman et al.lll980l ). With 
the shrinkage of the semi-major axis, the gap between the circumbinary disk and the binary 
would be reformed. The timescale for gap reformation can be estimated as r gap ~ 10" 4 r, 



4^cbd 

vis 
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using equation ([T]) (cf. lArtymowicz fc Lubowlll994l ). If r gr < r gap , mass transfer from the 



circumbinary disk will be stopped, and the light curves from the system would exhibit little 
periodicity as in Figure 4b. On the other hand, if r gr > r gap , mass will continue to be 
supplied from the circumbinary disk to each black hole, and the BBH system will maintain 
triple disks. In this case, the periodic light curves predicted for the accretion disks can serve 
as the electromagnetic counterparts to gravitational wave events detectable with instruments 
such as Laser Interferometer Space Antenna. 
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Fig. I. — Orbital phase dependence of the mass transfer rate M T from the circumbinary disk 
to the effective common gravitational radius of the black holes, and of the mass accretion 
rate M acc at the inner simulation boundary, r in = 5.0 x 10 _3 a, of the primary black hole. The 
right axis shows the bolometric luminosity Lboi corresponding to the mass transfer and mass 
accretion rate with an energy conversion efficiency rj = 0.1, normalized by the Eddington 
luminosity LEdd for a total black hole mass M B h = 1-0 x 10 8 M Q , where rj is defined by 
L bo i = r]M Bll c 2 . 



- 15 - 



1 1 

O.OIpc t=1Q 

+ 

o 

+ 


t=10.04 

o 


t= 1 0.1 

• 

+ 

+ 

o 


t=10.18 

• 

+ 

O 


t=10.5 

+ 

o o 

+ 


t=10.75 

o 

+ 

o 



Fig. 2. — Density maps of the two accretion disks around the supermassive BBH rotating 
with P or b ~ 9.4 yr and e = 0.5, for the period 10 < t < 11. Each panel shows on a logarithmic 
scale the surface density contours over a range of 5 orders of magnitude. The white crosses 
indicate the positions of mass input (see text). The supermassive BBH is rotating in a 
counterclockwise direction. The time is shown in each panel, and a length scale is shown in 
the top-left panel. The number of SPH particles during this period ranges from ~ 138, 000 
to 152,000. 
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Fig. 3. — (a) Radial distribution of the azimuthally averaged disk temperature at t — 10.0 
and at t — 10.18. (b) SEDs of the accretion disk around the black hole at t — 10.00 
and t = 10.18. The integrated luminosity, L v = d 2 uF u , is normalized by the Eddington 
luminosity L^d for a total black hole mass Mbh = 1-0 x 10 8 M Q ; F v is the flux at frequency 
v and d is an arbitrary distance from the observer. 
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Fig. A.— Orbital phase dependence of the UV (v ~ 1.0 x 10 16 Hz), optical (B band, 
v ~ 6.8 x 10 14 Hz), and NIR (K band, v ~ 1.36 x 10 13 Hz) light curves emitted from the 
disk. Panel (a) shows the light curves from the simulation with mass transfer, and panel (b) 
shows the light curves from the simulation with no mass transfer. 
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